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Abstract 

A Hamiltonian cycle of a graph is a closed path that visits every vertex 
once and only once. It serves as a model of a compact polymer on a lattice. 
I study the number of Hamiltonian cycles, or equivalently the entropy of a 
compact polymer, on various lattices that are not homogeneous but with a 
sublattice structure. Estimates for the number are obtained by two methods. 
One is the saddle point approximation for a field theoretic representation. 
The other is the numerical diagonalization of the transfer matrix of a fully 
packed loop model in the zero fugacity limit. In the latter method, several 
scaling exponents are also obtained. 
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1 Introduction 



A Hamiltonian cycle of a graph is a closed self-avoiding walk which visits 
every vertex once and only once. The number of all the Hamiltonian cycles 
of a graph G is denoted by 7i[G). Hamiltonian cycles have often been used to 
model polymers which fill the lattice completely 0. The quantity log7i(G) 
corresponds to the entropy of a polymer system on G in the compact phase. 
One can enrich the model by introducing a weight that depends on the shape 
of cycles in order to model more realistic polymers. For example, the polymer 
melting problem is studied by taking into account the bending energy ||. 
Another extended model relevant to the protein folding problem is proposed 
in refs. [|, 0. The quantity 7i{G) is also related to the partition function 
of zero-temperature O(n) model in the limit n — > 0. 

For homogeneous graphs (lattices) with N vertices, one expects that Ti{G) 
behaves as 



H{G) -> C{G)N 1 ~ 1 u N (N -> oo), (1) 
where the entropy per vertex logo; is defined by 

logo; = Km ^logft(G). (2) 

The connectivity constant u is supposed to be a bulk quantity which does 
not depend on boundary conditions, whereas the conformational exponent 7 
depends on such details of graphs ||. I assume that eq.(|l|) is not multiplied by 
a surface term pS N "\a < 1 proposed in refs.JT], |^, [| because in the present 
case Hamiltonian cycles fill the whole lattice which itself does not have a 
boundary. 

A field theory representation of 7i(G) for an arbitrary graph G is in- 
troduced by Orland, Itzykson, and de Dominicis [H| and has been used to 



study the extended models [0, |5j as well as the original one|]TT], [L2|]. For 
homogeneous graphs with the coordination number q, the saddle point ap- 
proximation yields 

H(G) ~ (u S p) N , usp:=^ (3) 

or uj ~ o; S p. Eq. (|D has been very successful. For the square lattice, 
= 4/e = 1.4715 • • • is quite near to a numerical estimate u ~ 1.473 • • • 



1 



14] , jlB] . It is found that the estimate wsp = 4/e is unaltered in the next 



leading order in the perturbation theory fl2 |. For the triangular lattice, the 
field theory predicts u>sp = 6/e = 2.20728 ■ ■ • while a numerical calculation 
suggests u) ~ 2.095 ■ ■ ■ |13| . An exact solution is available for the hexagonal 



lattice 0, [H . It implies a = 3 3 / 4 /2 = 1.13975 - ■ •, which is near to the 



estimate u^p = 3/e = 1.10364- • ■ . More examples are given in ref. [!1|. The 
saddle point equation is so good that one may speculate that the system is 
dominated by the saddle point configuration. It is desirable to understand 
why this approximation works so well. 

In this article, I systematically study 7i(G) for inhomogeneous graphs 
by the field theory. A graph is said to be homogeneous if its automorphism 
group acts on the set of vertices transitively. Roughly, a graph (lattice) is 
homogeneous if all the vertices are equivalent. Note that the estimate (|3]) 
is not applicable to inhomogeneous graphs. An example of inhomogeneous 
graph is the square-diagonal lattice shown in Fig. |l|. An earlier study in this 
direction can be found in ref. [ |TT|] . 

Inhomogeneous graphs are as important as homogeneous ones in physics 
because they model certain realistic materials well. Applying the saddle point 
method to them is expected to cast light on the nature of this extraordinary 
good approximation. 

Another achievement in this work is an accurate measurement of the 
connectivity constants u and scaling exponents by the numerical transfer 
matrix method. This is based on the mapping onto the fully packed loop 
model in the zero fugacity limit. The results should be compared with the 
field theoretic estimates. 

The organization of the paper is as follows. In section 0, the field theo- 
retic representation of 7i(G) is explained in order to make the presentation 
self-contained. I derive a formula for an estimate of 7i{G) by the saddle point 
method when a graph G has a sublattice structure which makes G inhomo- 
geneous in section 0. It is applied for a number of examples in section ^. In 
section [|, I prove some exact relations among 7~C{G) for several lattices. I 
explain the method and the result of numerical transfer matrix analysis in 
section || I summarize my results in section [7]. 
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2 Field theoretic representation 



Let G = (V, E) be a graph, where V and E stand for the sets of vertices and 
edges. One sets V = {n, r 2 , . . . , r^}, N := j^V . The number of Hamiltonian 
cycles, denoted by 7i(G), can be written as a 2iV-point function of a lattice 
field theory in a certain limit. To this end, one introduces O(n) lattice field 
(f>j(r) (r G V, j — 1, . . . ,n) with the action 



S[${r)] := i £ 0i(r)(V=TA 

r,r'6V,j,fc=l,... ,n 



-1 



(4) 



An infinitesimal parameter e — > +0 is the convergence factor. The N x N 
matrix A is the adjacency matrix of the graph G: 



1 if r, r' G V is connected by an e G £, 
otherwise. 



(5) 



It can readily be shown that JTU 



H(G) 



lim lim — 



2 (ri) 2 (r7v) 



(6) 



where 



{■■■)■, 



1 



n d ^ 

r£V,j=l,... ,n 



(7) 



with the normalization factor Z to ensure (1) = 1. 

The proof is diagrammatic. Note that the propagator is proportional 
to the connectivity matrix: {<j>j{r)(j)k(r')) = —y/—15jkA r r'- When Wick's 
theorem is applied, each of the surviving diagram after taking the n — > 
limit corresponds to a Hamiltonian cycle of G with an equal weight. The 
limit n — > is taken for discarding disconnected diagrams. 

When the graph G is homogeneous, there is a mean field saddle point 



r = 



(8) 



where q := ^ r A rr 
estimate @. 



is the coordination number of G. It gives rise to the 
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3 Saddle points for sublattices 



Given a graph G = (V, E) which is not necessarily homogeneous, I look for 
a saddle point for the integral expression (|). 

I decompose the set V into classes: V = U™ =1 V a , with iV a := #V a . Let 
g a (r) be the number of edges one of whose endpoints is r 6 V and the other 
belongs to V a . I only consider such divisions that q a (r) = q a (r f ) holds if r 
and r' belong to an identical class, say, V b . Then it makes sense to define 
q a b ■= q a (r). Because of the relation q ab N b = q ba N a , S ab := q a bN/N a becomes 
a symmetric matrix. 

Any graph does have such a decomposition; one is free to take V a = 
{r a }, a — 1, . . . , m = N. I concentrate, however, on the case where m stays 
finite in the limit iV — > oo. In other words, I focus on G with a sublattice 
structure^]. Hereafter V a is referred sublattice. 

I proceed to evaluating @ by the saddle point method. Respecting the 
sublattice structure, an ansatz for the saddle point configuration 0^(r) for 



(o) (r)=J>A(r) (9) 

a=l 

is employed. Here v a is an n-vector variable to be solved, labeled by a = 
1, . . . , m. The membership function 5 a (r) is defined by 

WH 1 ifrGK ' (io) 

otherwise. 

In order to calculate S*^^], one needs to know how A -1 acts on 5 a . First, I 
write the the number of edges connecting V a and r e V in two ways as 

m N 

J2 A rrMr')=J2w^Mr)- (11) 

r'eV 6=1 b 



1 Note, however, the saddle point method is worth considering as an approximation 
scheme even for the case V a — {r a }, m = N. Estimating H(G) by the saddle point method 
as described below takes only polynomial time in m(= N) while determining Ti(G) exactly 
for an arbitrary graph G is an NP-complete problem J17[ . 
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Rewriting the right hand side in terms of the m x m symmetric matrix S a b, 
I have 

TV m 

— &rr>W) = SaMr)- (12) 

a r'eV 6=1 

I apply A _1 S _1 on both sides to obtain A~ 1 5 a : 
m N 

E^U^r) = J2(A^) rr ,5 a (r')- (13) 

6=1 b r'GV 

Because of this relation, S[<p^(r)] reduces to 

S[fi°\r)] = -N\n2 + J2N a \n(v a ■ v a ) - - ]T {S~ x ) ah (v a ■ v b ). (14) 

a=l a, 6=1 

One can assume that the saddle point is of the form v a = (x a , 0,0,... , 0) 
and vary the action (|i~4"|) with respect to x a . One obtains m saddle point 
equations 

2N m 

-j^ = ^2{S~ 1 ) ab x a x b , (a = l,...,m) (15) 

6=1 

or equivalent ly 

m ^ 

J^x^qbaX' 1 = -. (16) 

6=1 

which are to be solved for This is a set of m quadratic 

equations with m unknowns. If the solution to ( TT5] ) is not unique, the one 
that minimizes the action fli~4]) should be selected. Putting the solution back 
into (PI), one obtains an estimate 

1 m / 2 \ N a /N 

H(G) ~ {u)sp) n , ^sp:=-n(y) • ( 17 ) 

a=l ^ ' 

This generalizes the estimate @ to the case of inhomogeneous graphs. 



5 



Figure 1: The square-diagonal (Union Jack) lattice. The decomposition of 
V consists of two classes: V\ = {•} and V 2 = {o}. 



4 Decorated square lattices 
4.1 Square-diagonal lattice 



I apply the estimate (|17D to the square-diagonal (Union Jack) lattice drawn 
in Fig. |l|. The decomposition of V consists of two classes V% = {•} and 
V 2 = {o}. The matrix q ab is given by 

'=(40)- < 18 > 

Each class has an equal number of vertices: N\ = N 2 = |iV. I assume that 
the sublattice structure is respected at the boundary. 

The field configuration x% = 8/5,x 2 = 5 becomes a saddle point in the 
limit 5 -> 0. Eq. yields 



= \ (19) 

which is, surprisingly, exactly the same value as the simple square lattice. 
It means that the T~C(G) of the square- diagonal lattice grows no faster than 
that of the simple square lattice though the former has many additional 
diagonal edges. It turns out that it is the case; the simple square lattice and 
the corresponding square-diagonal lattice have exactly the same number of 



Hamiltonian cycles, which I prove rigorously in subsection RTT . 
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Figure 2: Square-triangular lattices V(t = 2, it = 1) and D(t — l,u = 4). 
The integers t and u are the numbers of plaquettes with and without diagonal 
edges, respectively, consecutive in the horizontal direction. The sublattices 
are distinguished by the shape of the vertices. 



4.2 Square-triangular type lattices 

As another example of the use of (|I7]) , I apply it to a series of inhomogeneousf] 
square-triangular lattices D(t,u) and V(t,u) depicted in Fig. |^. They are 
square lattices with some of the diagonal edges added. The matrices q for 
them are all tridiagonal matrices if one labels sublattices appropriately. Here 
I show the m x m matrix q and the fraction N a /N of V(t = 1, u — 2m — 1) 
and D(t = 1, u = 2m - 3), m = (1), 2, 3, . . . . 

Lattice V(l,2m - 1) 



/4 1 
12 10 



12 
' ' • 



Vo 



••10 

' • . 2 1 

13, 



N a 1 

— = — (a = l,...,m). 

JM m 



(20) 



2 The lattice V(l, 1) turns out to be a homogeneous graph with the coordination number 
q = 5. 
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Lattice V(l,2m-2) 





/4 1 




\ 
• \ 




12 10 








1 2 






Q = 










o ' • • 


1 









2 


2 




V o 


1 


2/ 


Lattice V(2m - 1, 1) 








/ 3 2 




\ 
• \ 




2 2 2 








2 2 






9 = 










o ' • • 


2 









2 


4 




V o 


2 


2/ 


Lattice V(2m - 2, 1) 








/ 3 2 




\ 




2 2 2 








2 2 






Q = 










o ' • • 


2 









2 


2 




V o 


2 


4/ 


Lattice D(l,2m-3) 








/ 2 2 








4 2 








2 






9 = 










' ' • 


2 












4 




V o 


2 


O/ 



AT 



iV 



Na 

N 



Na 

N 



2^1 (a = m), 
2^T (otherwise). 



1 

m 



,m). 



2m-l 
2m-l 



(a = mj, 
(otherwise). 



2m-2 
2m-2 



l,m), 



(otherwise). 
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Lattice D(l,2m-2) 



/2 2 

4 2 

2 

' • • ' • • 



Vo 







••02 

2 2, 



Na 

N 



2m- 1 

_2 

2m-l 



(a = l), 
(otherwise). 



(25) 



Lattice L>(2m - 3, 1) 

/ 2 
4 2 2 



Vo- 



2 2 
' ' • 



••20 

' • • 2 4 

2 2, 



Na 

N 



2m-2 
_2 

2m-2 



(a = 1, m), 
(otherwise). 



(26) 



Lattice D(2m - 2, 1) 

/0 2 ■■■ 0\ 

4 2 2 : 

02 2 

' ' • ' • • 2 



Vo 



••22 

2 4 



N 



2m- 1 
_2 

2m-l 



(a = l), 
(otherwise) 



(27) 



For small m, eq. Ql5|) can be solved analytically. For example, V(l, 3) turns 
out to possess the saddle point 



^ = 5 = 4/2 + ^= 



2 3 
which implies 



v/3' 



(28) 



V(l,3) 
SP 



1 + 2v/3 



1.64225' 



(29) 



As m increases, one is forced to solve eq. fllS]) numerically. Fig. ^ illustrates 
the dependence of logc^sp of V(t,u) (o) and D(t,u) (o) on / := £/(£ + w). 
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0.9 
0.8 
0.7 
0.6 
0.5 
0.4 



V(t,u) (SP, field theory) 
D(t,u) (SP, field theory) 
V(t,u) (TM, numerical) 
D(t,u) (TM, numerical) 



0.3 







0.2 



0.4 



0.6 



/= 



t+u 



0.8 



Figure 3: The estimates for u are plotted. The saddle point approximation 
u SP for V(t,u) (o) and D(t,u) (o) with 2t + u < 16 are compared with 
the result of transfer matrix calculation c^tm ( D and x with an error bar, 
respectively). For details at / = 1/2, 1/3, see Table |I| 

Note that / = 0, 1 correspond to the simple square lattice and the triangular 
lattice, respectively. 

For the purpose of comparison, I also plot an estimate wtmj the result of 
numerical diagonalization of the transfer matrices for V(t, u) (□) and D(t, u) 
(x). The method to calculate cjtm is explained later in section |6|. Here I 
only mention that, under plausible assumptions, c^tm agrees with the true 
value u within the precision of the numerical work and the error originating 
from the finite size effects, though it takes considerable memory and CPU 
time to compute itQ. On the other hand, calculating usp is quite easy 

I argue that the estimates usp capture the qualitative feature of cjtm 
though there is a certain amount of difference in the numerical values. 

3 In the present analysis, the numerical diagonalization of sparse matrices of dimension 
< 10 6 has been carried out. 
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First of all, one notices that uj^p for Z}(1, 2m — 3) is saturated just like 
the square- diagonal lattice: 

w gi,am-3) = f ( m = 2 ,3,...)- (30) 

This suggests that the diagonal edges contribute almost nothing to H(D(1, 2m— 
3))- 

In accordance with this fact, the numerical estimate uj^' 2m ^ also stays 
identical with iutm of the simple square lattice. I will argue in subsection |5]^ 
that uj for D(l, 2m — 3) and that for the simple square lattice should in fact 
be identical. 

In the field theory, the saddle point for D(l, 2m — 3) is given by 

x 2a+ i = 8/5, x 2a = 5 ($->0). (31) 

It is intriguing that one needs the limiting procedure 5 — > whenever one 
obtains the saturated value ujsp = 4/e. 

Secondly, apart from the series Z}(1, 2m — 3), the estimate uj^p increases 
with / for both V(t,u) and D(t,u) as naturally expected. The result for 
wtm confirms this naive expectation. Moreover, in both analysis, ujsp/tm for 
D(t,u) is always slightly lower than V(t,u) for a fixed (t,u). 

Lastly, some detailed structures of c^tm are reproduced in u; SP . For exam- 
ple, there is a violation of simply-increasingness of ujtm for the pair .0(2, 1) 
and .0(3, 2). It is present also for ujsp. 

Near / = 1, there is a considerable difference between the numerical 
values of ujsp and ujtm- Note, however, this difference is already there for 
the triangular lattice, for which the mean field saddle point approximation 
@ has been applied. Thus, this difference does not imply the fault of the 
extension (p~7|) . 

5 Exact relations among TL(G) 

I prove some equalities among H(G) for the square-diagonal lattice and the 
square-triangular type lattices. 

5.1 Square-diagonal lattice 

Let Gsd be a square-diagonal lattice of a rectangular shape with L x x L y = N 
vertices. The periodic boundary condition is imposed across the edges of the 
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rectangle giving Gsd the toric topology. I assume that L x and L y are even 
and the square- diagonal lattice structure is consistent with the boundary 
condition. I define a simple square lattice Gss as the lattice obtained by 
removing all the diagonal edges from Gsd- It also has N vertices. I shall 
prove H(Gst>) = T~C(Gs$). In other words, I will show that Hamiltonian cycles 
on Gsd passes vertical and horizontal edges but not diagonal edges. 

I decompose the set of vertices of Gsd into V\ and V 2 as depicted in Fig. |l|. 
Then N\ = N 2 = N/2. On a Hamiltonian cycle, there are N vertices and ./V 
edges. Suppose edges out of iV are diagonal ones. The two ends of a 
diagonal edge belong to V\. On the other hand, there are N — vertical 
or horizontal edges each of which connects V% and V 2 . Thus I have 

N x = (2 x N w + 1 x (JV - iV (c °)) x - (32) 

and 

N 2 = (N-N®)xl (33) 

which imply = 0. 

Because H(Gst>) = H(Gss) f° r Gsd and Gss at any finite size, one can 
take the limit iV — > oo to have u) SB = u ss . 

5.2 Lattice D( 1,2m - 3) 

Let G m be a finite D(l,2m — 3) lattice of a rectangular shape. Again I 
assume the periodic boundary condition on both directions, respecting the 
D(l, 2m — 3) structure at the boundary. I shall show TC(G m ) = H(Gss)- 

For brevity, I write G C G' if V = V and E C E', where G = (V, E) and 
G' = {V, E'). Apparently, H{G) < H(G') if G C G'. We have seen that the 
inequality H(Gss) < ^(Csd) for the pair Gss E Gsd is actually an equality. 

It is crucial to find that Gss C G m C Gsd- Thus, I obtain H(Gss) = 
TC(G m ). This relation, in the limit N — > oo, again implies that a; ss = 

w £>(l,am-3)_ 

6 Numerical transfer matrix method 

In order to estimate u numerically, I map the problem of Hamiltonian cycles 
onto the zero fugacity limit of the fully packed loop model. Then I represent 
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the fully packed model onto a state sum model with a local weight in order to 
construct a row transfer matrix. The quantity u is related to an eigenvalue 
of the transfer matrix. 

The partition function of the fully packed model with the the loop fugacity 
n is given by 

Z FPL (n)= £ n*>. (34) 

fully packed loop configuration 

The sum in ([34]) is over all the fully packed loop configurations (Fig. |]), that 
is non-intersecting closed loops on the lattice such that every vertex is visited 
by exactly one of the loops. The number of connected components of loops, 
denoted by A/l, can be greater than one. 



Figure 4: An example of fully packed loop configuration with A/l = 8 on the 
simple square lattice V(0, 1) . The loops are drawn in thick lines. 



One can construct an equivalent state sum model in the following way. 
One begins with the triangular lattice V(l, 0) with the coordination number 
q = 6. The local degrees of freedom z live on each edge e. The three possible 
values of z(e) is <— (a directed edge), — > (an oppositely directed edge ), and 
— (a vacant edge). Let ei, . . . , e§ be the edges that share a vertex r. States 
on ei, . . . , e§ interact on r. The local vertex weight is 

W(z( ei ),... ,z(e 6 )) = s k (sgC) (35) 

if there is exactly one ingoing arrow and an outgoing one at r, where the 
integer k is given by 

(the angle of the right turn) 
k = . 36 

7T/4 
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Figure 5: A vertex weight for the state sum model for the fully packed 
loop model. The parameter s is a complex number. Other configurations are 
prohibited ( weight W = ). 

as illustrated in Fig. [5[ Otherwise, W = 0. 

In other lattices D(t,u) and V(t,u), some of the diagonal edges in the 
lattice V(l, 0) are missing. One extends the above construction to them by 
just interpreting these missing edges as existing edges on which only the 
vacant edge state is allowed. 

The partition function of the state sum model is given by 

z (e)=+-,-,-> (eG_E) reV 

Edff understood to share a vertex r. One sees that the sur- 

viving configurations in the state sum model are nothing but the fully packed 
loop configurations with a direction associated with each loop component. 
Let us inspect a contribution from a loop. When one walks along the loop in 
the associated direction to come back to the original point, one has changed 
the direction ±27r in total (if it is a topologically trivial loop). Therefore, the 
product of the vertex weight on the loop is s ±s due to the choice (^6|). This 
is why the partition function of the state sum model agrees with (|34]) with 

n = s 8 + s~ 8 . (38) 

Now one can think of a row transfer matrix T of the fully packed loop 
model. For the lattices D(t, u) and V(t, u), one introduces a transfer matrix 
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Figure 6: Square-triangular lattices V(t = 2,u = 1) and D(t = l,u = 2). 
The periodic boundary condition in the horizontal direction is introduced to 
give the lattice infinite cylinder topology. The seam is indicated by broken 
lines. 

which maps a state on a row of vertical and horizontal edges onto that on 
the upper row in Fig. 0. That is, a state acted by T is a linear combination 
of an arrow configuration on the row \/\/\/\/ ■ ■ ■ /\/\- For the lattice D(t, u), 
a one-unit shift in the horizontal direction is included in T to take care of 
the change of the positions of diagonal edges. 

It is important to note that the transfer matrix T commutes with the 
operator giving the net number of upward arrows 

d=(#T) + (#/)- (#!)-(#/)• (39) 

Thus T is block-diagonalized as T = ©^T^. 

One considers an infinite cylinder geometry which is suitable for numerical 
transfer matrix calculation. The horizontal direction is compactified with the 
period L as indicated in Fig. |6]. In the d = sector there can be loops which 
wind the cylinder once. It gives rise to a complication because such a loop 
contributes s° + s° = 2, not s 8 + s~ 8 = n. To avoid this, one introduces a 
seam as in Fig. || One declares that a loop which goes across the seam from 
the left to the right should gain an additional weight s 8 , while the left-going 
one should gain s~ 8 . 

Hamiltonian cycles are encoded in the d = sector in the limit n — > 0, 
or s -> exp[±Y|] on the infinite cylinder. The condition d = excludes the 
configurations which have unbalanced strings which travels from an end to 
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the other end of the cylinder. The limit n — > excludes small disconnected 
loops. In the following, n is simply set to zero. The connectivity constant 
ujtm{L) for this geometry is given by 

\ogLU TM (L) = j\og\\ (L)\, (40) 

where \ l d (L) is the i-th largest (in the absolute value) eigenvalue of T<f(L). 
One expects that, by taking the L — > oo limit, one arrives at the bulk value 
wtm = ^tm(°o)- The values of u TM in Fig. || have been calculated in this 
way. 

I notice that that the relation cjf^ = c^m = ™ 3 ^ holds already at 
every finite L. 

I assume that the rotational symmetry is restored in the L — > oo limit 
with the sound velocity v — 1 and that the system is described by a conformal 



field theory [jig] , which is the case for the simple square lattice V(0, 1). For 
other lattice, one observes, at least, the mass gap closes at L —>■ oo. Then 
scaling exponents are related to the finite size behavior of other eigenvalues 



of the transfer matrix ]T5| , pU| . The central charge c appears asj 

logu TM (L) = T log |A°(L)| = logWTM(oo) + — - + O(L^). (41) 

The correlation lengths & and the scaling dimensions Xi of general geometric 
operators labeled by i are given by [TJ] 



In particular, the geometric scaling exponents X± and X 2 (corresponding to 
one and two spanning strings, respectively) are related to the conformational 
exponent 7 in eq. ([!]) by 

7 = 2(1-X 1 )z/, i = 2-X 2 . (43) 

The result of the finite size scaling analysis is shown in Table [I]. I have 
treated even and odd L separately for V(0, 1) because there is an oscillatory 

4 For the triangular lattice, the restoration of conformal symmetry occurs in the frame 
where the triangles are regular. Thus the formulae should be modified by the factor of 
V3/4. 
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1 -L-L' 

lattice 




c 


Xi 


log w TM 


log^sp 


V(0, l)L=0mod2 


14 


-1.000(11) 


-0.04423(81) 


0.387171(1) 


0.386294 


V{0, l)L=lmod2 


13 


-2.510(5) 


-0.1875(10) 


0.387168(2) 


0.386294 


0)L=0mod3 


9 


-2.06(8) 


-0.196(18) 


0.74088(7) 


0.791761 


V(l,l) 


12 


-1.98(9) 


-0.175(17) 


0.563057(18) 


0.609438 


V(l,2) 


12 


-2.00(4) 


-0.189(3) 


0.498809(3) 


0.532186 


D(l,2) 


12 


-1.98(10) 


-0.186(25) 


0.470379(3) 


0.495166 


V(2, 1) 


9 


-1.93(9) 


-0.183(19) 


0.619823(4) 


0.665920 


D(2, 1) 


9 


-1.98(3) 


-0.189(14) 


0.571025(1) 


0.617343 


V(l,3) 


12 


-1.90(29) 


-0.176(3) 


0.46905(4) 


0.496068 


V(2,2) 


12 


-1.89(18) 


-0.216(10) 


0.55761(4) 


0.612848 


D(2,2) 


12 


-1.90(14) 


-0.203(10) 


0.53917(3) 


0.574521 


V(3, 1) 


8 


-1.87(18) 


-0.174(10) 


0.650(15) 


0.698199 


D(3, 1) 


8 


-2.01(18) 


-0.187(10) 


0.609(16) 


0.653694 


V(l,4) 


10 


-2.09(11) 


-0.187(4) 


0.452(11) 


0.474064 


D(l,4) 


10 


-2.10(11) 


-0.177(4) 


0.427(11) 


0.442671 


V(2,3) 


10 


-1.94(11) 


-0.170(3) 


0.522(10) 


0.553356 


D(2,3) 


10 


-2.08(11) 


-0.189(3) 


0.492(11) 


0.524924 


V(3,2) 


10 


-2.02(11) 


-0.189(3) 


0.595(11) 


0.634758 


D(3,2) 


10 


-1.99(11) 


-0.185(3) 


0.579(10) 


0.616846 


V(A, 1) 


10 


-1.95(11) 


-0.184(3) 


0.668(10) 


0.716732 


£>(4, 1) 


10 


-1.92(11) 


-0.180(3) 


0.636(10) 


0.683031 


\/(l,5) 


12 


-2.24(8) 


-0.213(3) 


0.441(8) 


0.459443 


1^(2,4) 


12 


-2.19(8) 


-0.210(3) 


0.499(8) 


0.525490 


D(2,4) 


12 


-2.07(8) 


-0.192(3) 


0.482(7) 


0.508479 


\/(3,3) 


12 


-2.04(8) 




0.559(8) 


0.593499 


D(3,3) 


12 


-2.00(8) 




0.530(8) 


0.564561 


^(1,6) 


14 


-2.02(6) 




0.432(5) 


0.459443 


D(l,6) 


14 


-2.10(6) 




0.411(6) 


0.386294 



Table 1: Exponents and the entropy per vertex estimated by the numerical 
transfer matrix method. Lattices D(l, 2m — 3) are suppressed because they 
share the properties with V(0, 1). The integer L max is the largest L, the 
size of the lattice in the horizontal direction, which has been examined. The 
exponent X 2 is also measured and is found to be zero for all the lattices. 



17 



behavior of period 2 due to the twist-like operator insertion |nj . I also see a 
period-3 oscillation for V(1,0). Presence of such oscillations is not assumed 
for analysis of other lattices. I find that the exponent X 2 is always zero at 
finite L and is supposed to be so in the limit L —>■ 00. 

The result in Table || suggests that, except for V(0, 1)(= D(0, 1)) and 
D(l,2m — 3), the system lies in the same universality class as the dense 
phase of O(n) loop model at n = whose exponents are 

c = -2, Xt = ~ = -0.1875, X 2 = 0, 

16 ! 19 (44) 

u = -, 7 = — = 1.1875. 
2 ' 16 

This fact can be understood well by regarding the fully packed loop model 
as a special case of the O(n) loop model whose partition function is given by 

Z loop (n,x- 1 )= £ x^- N n^, (45) 

loop configuration 

where the summation is over all the non-intersecting loop configurations not 
necessarily fully packed. Thus, the number Afy of vertices visited by a loop 
can be different from N. In the two parameter space (n, x _1 ), the dense phase 
emerges in the region where the O(n) temperature is small enough while 
the fully packed loop model is reproduced by setting x _1 to exactly zero. 

For the simple square lattice and the hexagonal lattice, the line x^ 1 = 
is an unstable critical line where a new universality emerges. This is due 
to the symmetry x^ 1 <-> — x^ 1 . Actually, this symmetry holds because any 
loops on the simple square lattice visit even number of vertices. 

The result in Table [I] means that x~ x = is not special for other lattices. 



This is because they admit loops which visit odd number of vertices [13 



7 Summary 

I have estimated the number Hamiltonian cycles on inhomogeneous graphs 
analytically and numerically. To estimate it analytically, I have employed 
the field theoretic representation and have applied the saddle point method. 
A formula for uj for graphs with sublattice structures has been obtained. The 
numerical estimation is based on the diagonalized the transfer matrix of the 
fully packed model on the infinite cylinder geometry. 
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The former result is simple and is believed to capture the physics of 
Hamiltonian cycles. The latter is accurate and provides with the scaling 
exponents though it spends lots of computer time. The results agree each 
other qualitatively. It is confirmed that the success of the estimate u>sp = 
3/e, 4/e and 6/e by the former method was not accidental; it now works fine 
for a family of lattices yielding various values of u>sp- Moreover, the former 
method successfully predicts the relation to ss = cu SD = ( X ;- D ( 1 ' 2m ~ 3 ) . 

The relation u^p < u^p holds for all the pairs G C G' I have examined 
so far. It may be proved that this holds generally. 

The latter method is useful in calculating the exact value of H.{G) of G 
with the planar or cylinder topology. There seems to be, however, no simple 
way to apply this method to the lattices with the torus and the higher- 
genus topology, and three-dimensional lattices. This is due to the presence 
of numerous topological sectors of self-avoiding loops on the lattice |2TJ . In 
contrast, the field theoretic approach was able to predict a boundary con- 
dition dependence of T~t(G) for a family of toric lattices [0. Therefore two 
approaches are considered complementary each other. 
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